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TECHNICAL MEMORANDUM 


PARAMETRIC STUDY IN WELD MISMATCH OF LONGITUDINALLY 
WELDED SSME HPFTP INLET 

SUMMARY 


Welded joints are an essential part of pressure vessels such as the space shuttle main engine 
(SSME) turbopumps. Defects produced in the welding process can be detrimental to weld perform- 
ance. Recently, review of the SSME high pressure fuel turbopump (HPFTP) titanium inlet x rays 
revealed several weld discrepancies such as penetrameter density issues, film processing dis- 
crepancies, weld width discrepancies, porosity, lack of fusion, and weld offsets. Currently, the 
sensitivity of welded structures to defects is of concern. From a fatigue standpoint, weld offset may 
have a serious effect since local yielding, in general, aggravates cyclic stress effects. Therefore, the 
weld offset issue is considered in this report. Using the finite element method and mathematical 
formulations, parametric studies were conducted to determine the influence of weld offsets and a 
variation of weld widths in longitudinally welded cylindrical structures with equal wall thicknesses 
on both sides of the joint. From the study, the finite element results and theoretical solutions are 
presented . 


INTRODUCTION 


As the application of welded titanium structures in the aerospace industry is increasing, 
extensive experimental data are required to characterize typical weld defects encountered in produc- 
tion and to determine their effects on weldments. Recently, review of SSME HPFTP titanium inlet 
x rays revealed several discrepancies such as penetrameter density issues, film processing dis- 
crepancies, weld width discrepancies, porosity, lack of fusion, and weld offsets. The fabrication 
problems have been associated with the local stress levels due to these defects between the 
components from which the SSME HPFTP is manufactured. As a result of issues which have 
arisen concerning the structural analysis applied to evaluate the SSME HFPTP titanium inlet welds, 
an in-house task has been initiated at NASA/Marshall Space Flight Center. The purpose of the task 
is to (1) verify the adequacy of the Rocketdyne porosity specification as applicable to the HPFTP 
titanium inlet, (2) validate the Rocketdyne rationale and fatigue analysis for porosity utilized for 
material reviews (MR’s), and (3) perform a porosity sensitivity study at the most critical weld loca- 
tions to understand the effect of pore sizes, clusters, location of pores in welds, wall thickness, and 
weld offset on the structural adequacy of the inlet. The structural integrity evaluation of a welded 
joint is not a simple process because of complications due to component geometry, defects, and the 
gradient of material properties across the weld metal, heat-affected zone, and parent metal which 
result in nonuniform deformations. Although the structural integrity evaluation must include all 
these factors, only the weld offset issue is considered in this report. If pressure vessels are made of 
welded elements, there arises the possibility of local stress concentrations due to offsets caused by 
weld mismatch. These may give rise to relatively high local stresses which may lead to early 



failure of the pressurized structures. In this study, a case has been considered under the 
assumptions that the cylinder has a very' large radius, and weld mismatch is localized and only 
extends lor a short distance along the circumferential direction of. the cylinder. With these 
assumptions, the flat plate specimens containing an offset and being loaded in tension are con- 
sidered. The finite element analysis and theoretical analysis were conducted to determine the influ- 
ence of weld offset and a variation of weld width in a longitudinally welded cylindrical structure 
with equal wall thicknesses on both sides of the joint. As a result, the stress concentration factors 
due to weld offset in a flat plate loaded in tension with various mismatches are presented in this 
report. 


FINITE ELEMENT ANALYSIS 


In the current analysis, plate idealization ot the longitudinal weld in the straight cylinder, 
assuming that the applied load is resisted predominantly by the inplane stressing of the shell, is 
used for analyzing weld offsets subjected to axial loading in tension with a force P as shown in 
figure 1. Since the membrane stress is very' much larger than the bending stress as t/2r « 1 (r 
represents the radius of cylinder), the flat plate model could be considered [Ij. Mismatch between 
components of flat plates can result in substantial local increases in stress. That is, if a butt weld 
in a flat plate structure under tension is mismatched, an offset is produced which causes bending 
stresses at the weld in addition to the tensile stresses. Thus, to estimate the effects of parameters 
considered in weld offsets, an analysis was conducted using a flat plate idealization of the joint 
geometry at the section of weld offset. ANSYS plane strain finite element models [2] were con- 
structed to perform a parametric study to determine the effect ol weld ollsets. The linite element 
model includes weld metal and parent metal regions. The representative finite element model is 
shown in figure 2. The mesh is more refined in the vicinity of a junction between weld metal and 
parent metal because of the local variation in stress expected in this region. Dimensions of the 
analysis models to obtain an offset stress concentration factor ( K o]f ) are shown in table 1. In this 
table L, e. l, and u represent weld width, weld offset, wall thickness, and Poisson s ratio, respec- 
tively. Practically, it is recognized that the case of more than 100-percent mismatch appears. How- 
ever, in this study mismatches up to 60 percent have been tested. It is only considered that a flush 
weld bead would be at such a seam to act as a transition section, and that there is no discontinuity 
of meridian slope, but an eccentricity between the middle surfaces on the two sides of the discon- 
tinuity is included. Stress concentrations resulting from reduced thickness or reentrant comers are 
not considered in determining the stress concentration factors. 


THEORETICAL ANALYSIS 


If a long cylinder with very large radius-to-thickness ratios and mismatch between shell 
elements of the same thickness and slope is considered, the plane strain assumption could be made 
using a flat plate with a mismatched welded joint loaded in tension. Sechler (3] used the same 
assumption and showed that if the small element ABCD in figure 3 is assumed to be rigid and to 
rotate as a solid body under the action of the two plate elements attached to it. the problem 
resolves itself into two parts as follows. 
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(1) Moment and force equilibrium of the offset element ABCD 

If the offset is e and the length of the weld width is L, then the equilibrium of the element 
is given by 

2PR cos (ct + P) — 25fi sin (ct+P)-2A/ = 0 , (1) 

or 

PR (cos p cos a — sin (3 sin a) — SR (sin a cos (3 + cos a sin (3) — M = 0 , (2) 

where P, S, and M stand for axial tensile load, end shear, and end moment, respectively. 

Here 

sin P = -A- and cos p = , 

where R is the diagonal of the small element ABCD. 

Assuming a is small, 

sin a — a and cos a — 1 

Then equation (1) becomes 

P(e — La) — S(L + ea) — 2M = 0 , (3) 

(2) A beam column problem for the two sheet elements AG and CE, consisting of a canti- 
lever beam with an end shear S, an axial tensile load P, and an end moment M 

The cantilever beam is considered as shown in figure 4. The basic equation is 

El = Sx - P{\\) - y) - M , (4) 

dx- 

where the direction of positive moment is clockwise. 
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This has a solution 


v = 0| sinh k.x + i‘2 cosh k.x— +\o + ~- , 


Using the boundary condition that y - (dy/d.x) = 0 at v = / and that ElUry/d.x 2 ) 
= y„ at v — 0. 


y(.v) = 


/ — + -^~ tanh x\ sinh k.x (cosh k.x — 1) — 

\ Pk cosh X P j P p 


+ JL (X - tanh X) - -jj- ( 1 + sinh X tanh X - cosh X) , 
P A" * 


= -JL (X- tanh (1 + sinh X tanh X- cosh X) , 

Pk 


dv , , 

- 7 - (*) 

d.x 


S + -M- tanh X ) cosh kx-^- sinh lex- — , 

P cosh X P / P F 


4 - 10 ) L- - A + ^g>anhX . 

</.v P 1 cosh X J P 


From figure 3 it is observed that 


= Rol sin (3 = Ret 


a = -f- ( v = 0) 
d.x 
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(3) Equation for stress concentration factor due to weld offset 


The following equations are derived from equations (3), (8), and (10), 


P(c — La) — S(L + ca) — 2M — 0 


( 13 ) 


a - tanh X) — 0 + s ' nh * tanh cosh ^ 
7 Pk P 


a 


5. 

P 


/ 1 — cosh \ \ 
\ cosh \ / 


+ _M_ lan h\ , 


(14) 


(15) 


which are to be solved for a, S. and M in terms of P. Since the most serious problems arise from 
hiszh stresses (P large) and. in most practical cases, / is relatively large, then \ = l (PIEI)' ~ will be 
sufficiently large ( >4) so that the following approximations are justified: 


sinh A. -sr cosh X >> 1 


(16) 


tanh X -sr I 


( 17 ) 


Equations (13). (14). and (15) then reduce to 


P(e - La) - S(L + ea) - 2M - 0 


( 18 ) 


aL _ 
2 


_S_ 

Pk 


(X-D- 


M 

P 


(19) 


a = 


5 , Mk 
P P 


( 20 ) 


Solving simultaneously. 


S 


Mk 


A + 2 
A + 2X-2 


( 21 ) 


2 Mk X — 2 

P A + 2X - 2 


( 22 ) 
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where 


Pe = A/(A+ 2) 



2Me 

E/ 


X - 2 

(A + 2X-2) 



A = Lk = L y/piEI 


(23) 


(24) 


Detailed study indicates that the last term in equation (23) is negligible, therefore 


M - /*', . 

A + 2 

(25) 

/V* 

A + 2X-2 ’ 

(26) 

2 AdX-2) 

(27) 

(A + 2)(A + 2A-2) 


The maximum bending stress due to M is 


= 


( + /-) 


M Hi 2 ) 

/• 1 /J2 


( + /-> 



(28) 


The direct stress is 


P P 

cro = T = T • 


(29) 


where A is the area. 


Thus, the maximum total stress is 


(7/ = <T() + O'/, = -j- + QM— 

1 r 


(30) 


The stress concentration factor due to weld offset is therefore given by 


K„ it = = 1 + 4^- = I + 3(W/) 


CT() 


Pi 


(Lit) p+ I 


(31) 
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where 


p = k(t/2) 

= -j- VpJei 
= _l / [ 

2 V (Xr 3 )/( 12(1 — v 2 )) 


= \J (3(1 -v 2 )ct 0 )/£ . (32) 

Sechler’s expression (equation (31)) accounts for the reduction of transverse deformation due 
to in-plane tensile loading, often called the stress stiffening or geometric stiffening effect. This is 
more apparent if equation (31) is rearranged as 


ATon- = \+3(e/t)(\/(Lp/t+\)) , 


(33) 


where the Lp/t quantity, when factored on the elt term, accounts for the stress stiffening. Note that 
Sechler uses the plate flexural rigidity (ErV 12(1 -v 2 )) in the expression for p (equation (32)), 
making K o{{ a function of Poisson’s ratio, and that by including the stress stiffening effect, his 
expression is nonlinear. In this study, a strain (cro/E) of 0.25 percent was used as a minimal strain 
for the HPFTP inlet longitudinal welds. A higher strain would make the Sechler curves in the 
results section of this paper conservative. 

If the stress stiffening effects are ignored or assumed negligible, equation (33) reduces to 


E off = 1 + 3 {elt) , 


(34) 


which is the standard, small defection, linear formulation commonly used for E off estimation of 
longitudinal weld mismatch. 


RESULTS AND CONCLUSIONS 


Analyses were performed for various ratios of weld width ( Lit ) with 20-, 40-, and 60- 
percent weld offsets. Initially, the analysis was performed varying offset for a constant weld width 
{Lit = 2.86). As can be seen in figure 5, the increase in stress is almost proportional to the 
amount of mismatch. That is, if the weld offset is increased, then the stress concentration factor 
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(AT ofr ) is also increased as expected. The results were compared to the conservatively used 
results [4], Generally, if K otf is calculated using axial stresses, the values are much higher than 
those calculated using effective stresses. Also, a comparison shows that the currently used values 
are the highest and that the values from a finite element model with 60-percent offset has a good 
correlation with the theoretical value. 

In figures 6, 7, and 8, it can be observed that if the weld width is increased then K olt is 
decreased. For L/t= 1.0, K otf calculated by an axial stress component from the finite element analy- 
sis is higher than currently used r . Consequently, it is believed that an excessive weld width 
does not adversely affect the structural capability of the joint. A wide weld bead would have a 
lower A" o) | than a narrow bead of the same amount of offset. The wide welds are therefore struc- 
turally acceptable. 

The variation of Poisson’s ratio was also tested as shown in figures 5 and 9 with the same 
conditions since a formulation, K 0 rr= I +3 [eh), for the longitudinal weld cases is currently used 
[4]. This implies that the weld stiffness effects could be negligible. However, as can be seen in 
table 2 and figures 5 and 9, if the effective stresses are considered, the values of cases with 0 = 0.3 
are very different from the models with v= 0.0. On the other hand, if axial stresses are compared, 
the differences are minor. In theoretical analysis, since the values are calculated based on the axial 
stresses, the theoretical values show the same characteristics as the cases of the finite element axial 
stresses. However, it is noted that the values from the effective stresses with v = 0.3 are lower than 
those with u = 0.0. But, the values from axial stresses and theoretical stresses with u = 0.3 are 
higher than those with v = 0.0. The comparison was also made for the different offset ratio (e/t) 
with a various ratio of weld width-to-wall thickness (Lit) for v = 0.0 and u = 0.3 in figures 6 to 8 
and 10 to 12. As can be observed in these figures, regardless of a change of Poisson’s ratio and an 
amount of offset, the results show that as the bead width is increased the stress concentration factor 
is decreased. It should be noted that if the weld width is increased, accuracy of the currently used 
results is decreased Consequently, since the maximum distortion energy theory (Von Mises yield 
criterion) to predict onset of plastic flow is the most accepted yield criterion for a material with a 
reasonable amount of ductility, the effect of Poisson’s ratio has to be taken into account in calcu- 
lating the stress concentration factors in terms of effective stresses for a longitudinal weld case. 

The results of the material discontinuity study between the weld metal (£j) and the parent 
metal (E 2 ) are presented in figures 13 through 16 with a ratio of E\/E 2 equal to 1.0. 1.1. 1.2, and 
1.3. For this study, the models were constructed with variations of Lit equal to 1,2. 3. and 4 for 
a constant weld offset (40 percent). As can be seen in these figures, the effect of a material differ- 
ence between the weld area ( E 2 ) and the parent metal area (£j) could be considered negligible. It is 
once again observed in these figures that for small L!t ratios the stress factors obtained by finite 
element analysis are more conservative than currently used values. 


DISCUSSION 

In order to get more accurate stresses in weldments by the finite element analysis, the con 
verged finite element solutions could be considered with increasing the number ot elements. How 
ever, since the models for this study have the same element aspect ratio and the stresses con- 
tributed by the mismatch are usually maximum at or very near the joint, the information obtained 
herein is believed to be all that is necessary for the parametric study purposes. Nevertheless, the 
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nonlinear analysis including both geometric and material nonlinearities might be necessary because 
the potentially critical locations with respect to the local stress raisers, for example, through sharp 
notches, weld metal porosity, slag or oxide inclusions, poor penetrations, shrinkage, weld surface 
roughness, etc., are expected to occur at the weld-parent metal interface. Also, if a change in 
curvature occurs at the mismatch, additional discontinuity stresses result. In this case, the mismatch 
stresses should be superimposed upon the geometrical discontinuity stresses from the curvature 
change and membrane stresses to obtain the total stress. In addition to these stresses from the 
geometry changes, since a weld inherently produces metallurgical notches, this metallurgical notch 
due mostly to the heat effect in the heat-affected zone also needs to be considered. Therefore, if 
the precise assessment for the specific design problems is needed in weld discrepancies, all of these 
factors might be necessary as a whole in the analysis. However, generally speaking, for the 
parametric study to see an effect of the influence of a specific parameter the present analysis 
models are considered to be fairly reasonable. 

It should again be mentioned here that the present analysis results are applicable only to 
shells of constant thickness with a large ratio of radius-to-wall thickness, and stresses are found 
only at the junction. Finally, the information provided in this study should be of help in the selec- 
tion of allowable amounts of mismatch for improving design and manufacture of the SSME 
HPFTP. 



Original 

Deflected 

(a) Detlection of middle surface 



(b) Loading and geometry 

Figure 1. Weld offset of longitudinally welded cylindrical structures [5], 



Figure 2. Representative finite element model of weld offset. 
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Figure 3. Forces 



Figure 4. As; 
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Model # L e t L/t e/t U E1 / E2 

1 0.2 0,014 0.07 2.86 0.2 0 3 1.0 

2 0.2 0.028 0.07 2.86 0.4 0.3 1.0 

3 0.2 0,042 0.07 2.86 0.6 OJ 1.0 

4 0,07 0.014 0.07 1.0 02 0.3 1.0 

5 0.14 0.014 0.07 2.0 0.2 0.3 1.0 

6 0,21 0.014 0,07 3.0 0.2 0.3 1.0 

7 0.28 0.014 0.07 4,0 0,2 0.3 1.0 

8 0.07 0.028 0.07 1.0 OA 0.3 1.0 

9 0.14 0.028 0.07 2.0 04 0.3 1.0 

10 0,21 0.028 0.07 3.0 04 OJ 1.0 

11 0,28 0.028 0.07 4,0 0.4 0.3 1.0 

12 0.07 0.042 0.07 1.0 0,6 03 1.0 

13 0.14 0.042 0.07 2.0 0.6 0.3 1.0 

14 0,21 0.042 0.07 3.0 0.6 0.3 L0 _ 

15 0.28 0.042 0.07 4J ) 0.6 0.3 1.0 


0.2 0.014 0.07 2.86 I 0.2 I 0.0 I LO 


17 

0.2 

0.028 

0.07 

2.86 

0.4 

0.0 

1.0 

18 

0.2 

0.042 

0.07 

2.86 

0.6 

0.0 

1.0 

19 

0.07 

0.014 

0.07 

1.0 

0.2 

0.0 

1.0 

20 

0.14 

0.014 

0.07 

2.0 

0.2 

o.o 

1.0 

21 

0.21 

0.014 

0.07 

3.0 

0.2 

o.o 

1.0 

22 

0.28 

0.014 

0.07 

4.0 

0.2 

■031 1 

1.0 

23 

0.07 


0.07 

1.0 

0.4 

0.0 

1.0 

24 

0.14 

0.028 

0.07 

2.0 

0.4 

0.0 

1.0 

25 

0.21 

0.028 

0.07 

3.0 

0.4 

0.0 

1.0 

26 

0.28 

0.028 

0.07 

4.0 

0.4 

0.0 

1.0 

27 

0.07 

0.042 

0.07 

1.0 

0.6 

0.0 

1.0 

28 

0.14 

0.042 

0.07 

■bob 

0.6 

0.0 

1.0 

29 

0.21 

0.042 

0.07 

3.0 

0.6 

0.0 

1.0 

30 

0.28 


0.07 

4.0 

0.6 

0.0 

1.0 

31 

0.07 

0.028 

0.07 

1.0 

0.4 

0.3 

1.0 

32 

6.07 

0.028 

0.07 

1.0 

0.4 

0.3 

1.1 

33 

0.07 

0.028 

0.07 

1.0 

0.4 

0.3 

1.2 

34 

0.07 

6.028 

0.07 

1.0 

0.4 

0.3 

1.3 

35 

0.14 


0.07 

2.0 

0.4 

0.3 

1.0 

36 

0.14 

0.028 

0.07 

2.0 

0.4 

0.3 

1.1 

37 

0.14 

! 0.028 

0.07 

2.0 

0.4 

0.3 

1.2 

38 

0.14 

0.028 

0.07 

2.0 

0.4 

0.3 

1.3 

39 

0.21 

0.028 

0.07 

3.0 

0.4 

0.3 

1.0 

40 

0.21 

i 0.028 

0.07 

3.0 

0.4 

6.3 

1.1 

41 

0.21 

0.028 

0.07 

3.0 

0.4 

0.3 

1.2 

42 

0.21 

0.028 

0.07 

3.0 

6.4 

0.3 

1.3 

43 

0.28 

0.028 

0.07 

4.0 

0.4 

0.3 

1.0 

44 

0.28 


0.07 

4.0 

0.4 

0.3 

1.1 

45 

0.28 

0.028 

0.07 

4.0 

0.4 

0.3 

' 1.2 

46 

0.28 

0.028 

0.07 

4.6 

6.4 

0.3 

1.3 





















































































































































































































































































































































Model 


1 


2 


3 


FEM(eff 


1.370 


1.886 


2.439 


1.483 


1.395 


1.358 


1.318 


2.009 


1.920 


1.833 


1.760 


2.535 


2.420 


2.315 


2.225 


1.535 


2.109 


2.724 


1.658 


.583 


1.521 


1.477 


2.222 


2.145 


2.044 


1.967 


2.806 


2.710 


2.580 


2.482 


2.009 


2.015 


2.022 


1.943 


1.948 


1.954 


1.961 


1.833 


1.835 


1.838 


1.842 


1.760 


1.764 


1.768 


1.773 


FEM(axial 


1.543 


2.128 


2.750 


1.677 


1.571 


1.530 


1.482 


2.239 


2.150 


2.058 


1.978 


2.803 


2.700 


2.598 


2.497 


1.539 


2.117 


2.731 


1.666 


1.586 


1.526 


1.479 


2.214 


2.150 


2.050 


1.972 


2.785 


2.690 


2.582 


2.486 


2.239 


2.251 


2.264 


2.277 


2.185 


2.193 


2.202 


2.210 


2.058 


2.062 


2.067 


2.074 


1.978 


1.981 


1.985 


1.989 


1.485 


1.970 


2.455 


1.544 


1.515 


1.480 


1.451 


2.108 


2.029 


1.961 


1.901 


2.662 


2.544 


2.441 


2.352 


1.481 


1.961 


2 


1.552 


1.511 


1.476 


1.445 


2.104 


2.022 


1.952 


1.890 


2.656 


2.533 


2.428 


2.366 


2.108 


2.112 


2 . 


2.119 


2.029 


2.036 


2.042 


2.048 


1.961 


1.970 


1.978 


1.985 


1.901 


2 


1.921 


1.930 


1+3 (e/t 


1.600 


2.200 


2.800 


1.600 


1.600 


1.600 


1.600 


2.200 


2.200 


2.200 


2.200 


2.800 


2.800 


2.800 


2.800 


1.600 


2.200 


2.800 


1.600 


1.600 


1.600 


1.600 


2.200 


2.200 


2.200 


2.200 


2.800 


2.800 


2.800 


2.800 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 


2.200 
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